Conductance Properties of Carbon-Based Molecular Junctions 
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We present a comprehensive study of the properties of the off-resonant conductance spectrum in 
oligomer nanojunctions between graphitic electrodes. By employing first-principle-based methods 
and the Landauer approach of quantum transport, we identify how the electronic structure of the 
molecular junction components is reflected in electron transport across such systems. For virtually 
all energies within the conduction gap of the corresponding idealised polymer chain, we show that: 

a) the inverse decay length of the tunnelling conductance is intrinsically defined by the complex- 
band structure of the molecular wire despite ultrashort oligomer lengths of few monomer units, and 

b) the contact conductance crucially depends on both the local density of states on the metal side 
and the realised interfacial contact. 
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I. INTRODUCTION 

Molecular electronics has been an active research field 
d H El since the early seventies, when a rectifier based 
on a single organic molecule was proposed |4| . Exten- 
sive experimental H H H H 
El B E , E|l and theoretical 

m m m mml m m m im 1351 m m m, \m 

investigations have examined the electronic response of 
a single molecule hybridising with metals which act as 
donor and acceptor reservoirs. Electron transport across 
metal/molecule/metal junctions (MMMs) depends cru- 
cially on the details of the contacts and the exact choice 
of materials, so that the fundamental mechanisms have 
been highly controversial. Fine tuning may lead to a 
range of behaviours in the current- voltage characteris- 
tics 0. Examples include negative differential resistance, 
strong rectification, molecular memory phenomena, and 
even the Coulomb blockade effect, a hallmark of meso- 
scopic physics. 

For a class of MMMs increasing consensus based on ex- 
perimental findings m El El El El El 

fli^ l points to 

a dominant transport mechanism, namely, to through- 
bond tunnelling. Such setups typically consist of rela- 
tively short oligomers like alkane [l3l IbH fTSL ITr3| and 
phenyl-based chains [TtI ITsl ITflj which are covalently 
bonded either directly or via anchor groups to the metal- 
lic electrodes. The metal of choice is usually Au for its 
noble properties. Most recently, junctions with graphitic 
and mercury drop contacts have also been investigated 


Experimental results show temperature independent 
electron transport and a conductance g which decreases 
exponentially versus the molecular length L, with an in- 
verse decay length (3 deriving from the oligomer. Namely, 



g(E F ,V) = go (E F ,V)e 



-f3(E F ,V)L 



(1) 



where Ep denotes the equilibrium Fermi energy. The 
pre-exponential factor g D is determined by the actual in- 
terfacial contact. In most realisations, measured quanti- 
ties do not depend on the source-drain voltage V for low 
bias. In another experiment |15|. a weak dependence has 
been indicated. Nevertheless, such behaviour is consis- 
tent with the picture of coherent off-resonant tunnelling 
across the molecular bridge. 

Eq. Q]has been theoretically investigated deep inside 
the tunnelling regime, both analytically for 7r-electron 
models of increasing complexit y 14 ll I42L l43l and numeri- 
cally for some realistic MMMs [29l I44L l45l| .In those stud- 
ies, electron-vibrational coupling has been disregarded 
and the zero temperature T, zero bias limit have been 
taken. Up to T not much lower than room tempera- 
ture (< 25meV), all assumptions constitute a reasonable 
approximation for specific molecular junctions. This de- 
pends on the molecular vibrational frequencies and con- 
formational barriers, electron resident times, and the 
non-equilibrium electrostatic potential profile. Moreover, 
these assumptions are in agreement with the above ex- 
periments, allowing for an electrical response that is pre- 
dominantly a property of equilibrium conditio ns, such as 
the band line-up or charge relaxation [2I |2j lH |51 0, 
IH|2(|, the optimised atomic configuration of the molecu- 



lar bridge |2 

Ella 



and the contact microstructure 



Electron transport across a molecular junction has a 
lot in common with the intramolecular or donor-bridge- 
acceptor (DA) electron transfer. The exponential de- 
pendence of the conductance in Eq. ^ is reminiscent of 
the well-studied McConnel's superexchange mechanism 
in the weak coupling limit |46l I47L |48| , where 



^| HrM | 2 F~ e -^. 
n 



(2) 
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A connection between the electron transfer rate hp a and 
the conductance has been given by Nitzan 49] in certain 
limits. Hda is the coupling between donor and acceptor 
electronic states at a potential surface crossing energy 
E, and F is the Franck-Condon factor. However, the 
continuum of states of the reservoirs and their constraints 



redefine the problem and have to be explicitly considered 
for a full account of the conduction properties. From 
this viewpoint, a transport calculation across the whole 
system is required. 

Mujica et al ^lj generalised previous analytical results 
of Everson and Karplus 0] on electron transfer across a 
single-orbital tight-binding linear homogeneous chain, to 
study conductance properties. Their studies revealed the 
exponential length decrease for energies far-off the energy 
spectrum of the molecular bridge. Onipko et al ^3] found 
an approximate solution for conjugated oligomers in the 
7r-electron model. By disregarding the effect of the elec- 
trodes to the molecular electronic structure in the tunnel 
limit /3(E) L ^> 1, their results closely resemble the elec- 
tron transfer problem. However, they provide a simple 
expression for g Q . Using extended Hiickel calculations, 
Joachim with Vinuesa Q and Magoga jl^ have studied 
numerically the conductance at Ep of several relatively 
long conjugated molecules, where Ep falls in the middle 
of the wire conduction gap. They established that the 
damping factor is a characteristic property of the wire 
and that g Q is related to the contact realisation of the 
bridge, in agreement with Onipko et al 0- Moreover, 
an analytic expression for (3 was derived within a simple 
model [43j and its properties were studied within a level- 
repulsion approach of Random Matrix Theory HH1 ■ 

Most recently, Joachim and Magoga 1 an d Tomfohr 
and Sankey I~>.jI | , provided independently approximate ex- 
pressions for the effective mass of electrons in MMMs 
tunnelling through the oligomer-induced barrier. They 
considered the ideal polymer limit and approximated (3 
to the spatial decay parameter of the wavefunction. This 
possibility has been conjectured in many of the afore- 
mentioned studies. Tomfohr and Sankey studied the lo- 
cal DOS decay parameter of octanedithiols between Au 
(111) surfaces and found good qualitative agreement with 
the imaginary part of the wavevector of wavefunctions in 
the forbidden energy domain of the corresponding one- 
dimensional crystal, namely, the complex-band structure 

mm 

In this paper, we examine off-resonant electron 
transport across molecular junctions in view of the 
new experiments and advances of theoretical con- 
cepts. In particular, we provide a comprehen- 
sive analysis of the entire conductance spectrum 
inside the Highest-Occupied-Molecular-Orbital-Lowest- 
Unoccupied-Molecular-Orbital (HOMO-LUMO) gap Eg 1 
of phenyl-based oligomers, as a function of the num- 
ber of monomer units. We focus on studies of oligo- 
phenyl-ethynyl (OPE), planar and non-planar oligo-para- 
phenyl (OPP) molecular wires covalently bonded to elec- 
trodes consisting of graphitic ribbons. First-principlc- 
based methods are employed to calculate the electronic 
structure of the molecular junction components, and the 
conductance in the Landauer approach of quantum trans- 
port. 

We find that the conductance is almost always an ex- 
ponentially decreasing function of length and not only 



for /3(E) L ^> 1. We observe tunnelling characteristics 
for ultrashort oligomer molecular junctions and/or en- 
ergies not far-off the edges of the conduction gap E g of 
the corresponding polymer (PPE and PPP, respectively) . 
The contact conductance can be approximately written 
as go(E) = rj^Tfjf(E), where T L /u describes the atomic 
coupling at the molecule/metal interface, and f(E) re- 
flects both the electrode and molecular spectral proper- 
ties. The inverse decay length is determined by the en- 
ergy spectrum of the polymer when extended to include 
complex wavevectors. A remarkable one-to-one corre- 
spondence to the imaginary part of the latter is found 
for all cases. A brief account of the complex wavevector 
mapping applied to PPE has been previously given by 
the authors Mathematically, the problem of scat- 

tering through finite periodic systems has been studied 
to some extent in Ref. |59l |. 

Since interfacial barriers in our systems are generally 
lower than the ones realised in the experiments, these 
results explain qualitatively the observations in molec- 
ular junctions of a handful of monomers. In addition, 
they acquire special importance when recognising that 
Ep may lie anywhere inside the gap. It is in fact com- 
mon knowledge that it is difficult to a priori locate Ep 
which crucially depends on the combination of anchor 
groups and electrodes. 

The choice of OPP and OPE as molecular bridges in 
our MMMs is dictated by the usual experimental se- 
tups. Graphene-like macromolecules have been synthe- 
sised |(3(j and studied theoretically (3^ as the active com- 
ponents of molecular junctions. Moreover, driven by 
a vision of purely carbon electronics experiments with 
pyrolitic graphite as one of the metals have been re- 
ported 18], for which the structures we study constitute 
a first approximation. On the other hand, the graphitic 
ribbons we consider have some special electronic states 
HH EE EE 03 Ell with large local density of states 
around Ep , which are responsible for additional features 
in the conductance. In analogy to other carbon-based 
systems |36j| in which they give much richer spectra, these 
states provide an extra channel for resonant transport 
within the HOMO-LUMO gap of the molecular bridge 
[6(| . As a function of the length of oligomers, they also 
compete against the exponential law Eq. ^ 

The structure of the paper is as follows. In Sec. ITU we 
briefly introduce the theoretical methods that we used, 
and the details of the investigated systems. In Sec. IIIII 
we present the electronic structure of the studied molec- 
ular junctions with emphasis on the properties relevant 
to the aimed transport. We also analyse the concept of 
complex-band structure. The conductance properties are 
discussed in Sec. IIVI Finally, we make some concluding 
remarks. 



II. COMPUTATIONAL FRAMEWORK 



A. Electronic Structure and Transport Properties 



for a wider bias range for voltages dropping mostly at the 
contacts. Finally, the transmission function is given by 
the Green's function method of quantum transport via 



The electronic structure of the oligomer as well as that 
of the electrodes is treated within an approach based on 
Density Functional Theory (DFT) in the local density 
approximation (LDA). By employing a linear combina- 
tion of atomic orbitals, the method falls into the class of 
tight-binding (TB)-DFT. It has been successfully applied 
to the calculation of properties of a range of materials 
including semiconduct ors, carbon nanotubes, fullerenes, 
DNA and proteins [H El El El E3 

The single-particle electronic Kohn-Sham eigenstates 
tpi of the system are expanded in a non-orthogonal basis 
set <f^(r — R^) taken as a valence basis localised at the 
ionic positions R M , namely 



T(E F ) = Tr 



r i G M r 2Gj\/ 



(6) 



!c> M (r-i? M ). 



(3) 



When this Ansatz is substituted in the Kohn-Sham 
equations for ipi , it yields a set of algebraic equations 



(4) 



The crucial step in deriving Eq. 0] is a further ap- 
proximation of the full many-body problem by a tight- 
binding Hamiltonian which is a function of atomic den- 
sities 68]. The a priori parametrisation scheme for the 
two-centre, distance dependent Hamiltonian, -ffJ^P, and 
overlap, ST„, matrix elements provides an efficient algo- 
rithm for electronic structure calculations. Explicit com- 
putational details are given in Ref . [67j ■ Even though the 
form of Eq. 0]is evidently that of the extended Huckel, 
all necessary matrix elements are determined without in- 
troducing empirical parameters. 

To calculate the conductance, we adopt the Landauer 
picture of quantum transport |56l l57| in combination 
with the TB-DFT. In this formalism, the current through 
a structure is expressed in terms of the transmission func- 
tion T(E, V). The latter denotes the probability of elec- 
trons with energy E to propagate coherently across the 
system between the two electron reservoirs. Reading ex- 
plicitly, 

2e f +co 

I=^J dE T(E,V)[f(E-m,T)-f{E-^,T)], 

(5) 

where f(E — fi, T) is the Fermi function, and /i 1)2 refer to 
the chemical potentials of the two electrodes (1,2). The 
factor 2 in the equality implies spin degeneracy. 

Here, we take the low bias limit and drop the voltage 
dependence of T, which in general is a function of the ex- 
act electrostatic profile. This approximation yields accu- 
rate results for V ~ O.leV j^. It is also expected to hold 



where G r M (E) = (ESjf - H^f - S)" 1 is the retarded 
molecular Green function at equilibrium. 

Using the Lowdin projection technique, G r M is dressed 
by a self-energy interaction X = Si + X 2 that takes into 
account the hybridisation of the molecular bridge with 
the metals. The spectral widths iTi 2 act as an effective 
coupling to the electronic states of the electrodes. £1,2 
and T12 are given by 



Ei, 2 (25) = 3{ 2 _ M (E) G\ >2 {E) J x , 2 _m(£) 



and 



T 1>2 (E)=i Ei, 2 (£) 



Ej j2 (i5) 



(7) 



(8) 



respectively. J1.2-M = ESi, 2 -m - Hf f_ M is the gen- 
eralised expression for a non-orthogonal basis, instead 
of the simple occurrence of the coupling matrix H^_ M 
between the electrodes and the molecule. The retarded 
Green function G\ 2 °f the leads is calculated by the San- 
vito et al algorithm [jl|, also extended to include the 
overlap between atomic orbitals. The above procedure 
mainly involves the manipulation of Nm x Nm matrices, 
where Nm is the number of atomic orbitals taken at po- 
sitions on the molecule. Exception is Eq. [7|that includes 
electrode degrees of freedom. However, interactions are 
of short-range due to the tight-binding approximation, 
and matrix dimensions never exceed ~ Nm- 

In what follows, we focus on generic properties of elec- 
tron transport across carbon-based molecular junctions. 
To this end, we study the conductance spectrum which 
is defined via 



S(E) 



2e 2 



-T(E). 



(9) 



It follows from Eq. within the linear response that the 
above expression at Ep strictly coincides with the con- 
ductance for very low bias. For finite but small voltages 
such that non-equilibrium effects are insignificant gener- 
alised formulae involving g(E) may be derived (cf |6(]). 
However, a complete investigation of the voltage depen- 
dence is out of the scope of the paper and is left for the 
future. We note the main effect of moderate-to-room 
temperature is to broaden features of the conductance at 
the ~ O.leV scale. 



B. System Specifics 

We consider all-carbon molecular junctions such as 
those depicted in Fig. ^ The electrodes are (zigzag) 



Transverse 




FIG. 1: A representative molecular junction among those 
studied is the structure of diphenyl-ethynylene (N = 2) 
bonded to two graphitic electrodes along either the zigzag 
(left) or the armchair (right) edge. Open circles denote hy- 
drogen atoms. See Sec. Ill Bl for more details. 



graphitic ribbons of the indicated width and infinite ex- 
tension along the carbon boundary of zigzag shape. By 
taking into account the weak inter-plane interactions in 
graphite, we isolate a single plane for simplicity. The 
molecule is a conjugated molecular wire of a few re- 
peated monomer units, N, also taken as our length unit. 
To demonstrate topological effects on the conductance, 
we initially examine two different types of junctions; one 
with the oligomer being attached to the electrodes along 
the zigzag edge ('transverse' configuration), and one with 
the wire connected to the armchair side of the ribbons 
('parallel' configuration), as shown in Fig. All dan- 
gling bonds at edge sites are saturated by hydrogen atoms 
(open circles in Fig. [IJ, so tnat tne Y do not contribute 
to the density of states around the Fermi energy Ep. 

The molecular bridges under investigation are formed 
by two chemically different oligomers, and conformation- 
ally varying: a) oligo-phenyl-ethynylenc, and b) planar 
and non-planar oligo-para-phcnylcnc which differs from 
OPE in that it does not have the triple carbon-carbon 
bond between the benzene rings. This makes OPP more 
flexible with energetically favourable conformation the 
non-planar one. Tetrahedral angles between benzene 
rings may var y fr om 45° for very small oligomers to 27° 
for polymers |73l l74j . At room temperature, ring rota- 
tional motion occurs between the planar and non-planar 
conformation. When the temperature is cooled down 
both possibilities may be realised in experiments with 
self- assembled monolayers [l^. Therefore, we consider 
several geometries as summarised in Fig. [5] 




FIG. 2: Configurations of oligo-para-phenyl molecular bridges 
studied in this paper. They are distinguished between odd 
and even TV, e.g. for (a)-(d), N = 5, and for (e), N — 6, re- 
spectively. The graphene ribbon electrodes contacted at the 
end-sites on both sides are in-plane with the page. The tetra- 
hedral angle for non-planar OPP ((c)-(e)) and the contact 
geometry (b) is 45°. 

The structures (molecular bridge+electrodes) are opti- 
mised in all calculations. To simulate the bulk electrodes, 
a constrained atomic optimisation of junctions such as 
those of Fig. ^is performed. Graphene atoms away from 
the oligomer/ribbon interfacial bond by approximately 
3 x hexagonal lattice constant remain immovable in the 
positions they assume in perfect graphene. The inner car- 
bon atoms, which look unsaturated in Fig. ^ are bonded 
to H to avoid spurious effects during optimisation. As in- 
put, initial single C-C bond lengths equal 1.39A, triple 
C-C bond lengths equal 1.2A, and C-H distances equal 
to 1.09A. All angles are at 120°. 

III. ELECTRONIC STRUCTURE 
A. Electrodes: Graphene Ribbons 

There has been an extensive interest in the properties 
of nano- graphitic materials recently. Specifically, ribbon- 
shaped graphite as in Fig. ^ is intensively investigated 
and the scope of this section is limited to pointing out the 
electronic properties that are crucial for the conductance 
of the systems we study. These are related to peaks in 
the density of states (DOS). Such features arise from the 
low dimensionality of the ribbons, the topology of their 
edges, and the hexagonal lattice. 

In Fig. Ola) we plot the density of states of a ribbon 
with width M = 4. The latter is given by the maxi- 
mum number of honeycomb cells that can be found in 
a cross-section. The ribbon DOS is contrasted with the 
DOS of bulk graphene (dotted line). It exhibits the char- 
acteristic van-Hove singularities of quasi-one-dimensional 




simplicity of the 7r-electron Hamiltonian. In those, the 
topological origin of the peak and its stability against 
disorder were also established. 

Finally, we note that for applications in molecular elec- 
tronics, one should, more precisely, look at the local DOS 
since the local electronic structure is most important. 
The overall DOS may be misleading. For example, elec- 
tronic states in a quasi-one-dimcnsional square lattice or 
simply those of the Anderson-Newns model of chemisorp- 
tion possess van-Hove singularities but exhibit a smooth 
local DOS. This owes to a cancellation of border zone 
anomalies (group velocity v g oc sin(ka)) by the longitu- 
dinal component of the wavefunction [3J|. Nevertheless, 
this does not apply to our discussion. Using the analogy 
of graphene ribbons to carbon nanotubes, and in agree- 
ment with the topological arguments, we employ a pre- 
vious result [34| stating that for states with zero group 
velocity outside the border zone, as in Fig. no such 
cancellation occurs. Indeed, we have calculated the local 
density of states at the carbon atoms on the electrode 
sides which are bonded to the oligomers for all studied 
molecular junctions. Unless stated, we have found small 
deviations from the DOS of Fig. |3| (see also inset of Fig. 
[5|, mainly caused by the oligomer chemisorption. 



B. Molecular Wires: Poly-Phenyl-Ethynylene 
(PPE) and Poly-Para-Phenylene (PPP) 



FIG. 3: (a) Density of states (DOS) for graphitic systems; 
dashed and solid lines correspond to graphene and a zigzag 
ribbon, respectively. The inset shows the DOS of the ribbon 
at energies around Ef. (b) Band structure of a zigzag ribbon 
of width M = 4. 



systems. Apart from that one around Ep, which we dis- 
cuss separately, their peak positions depend on the width, 
and they are present independent of the topology of the 
edges. 

The double peak of the ribbon DOS around Ep (in- 
set of Fig. [TJ isgeneric to graphene edges with a zigzag 
geometry |6lL |62|. For zigzag ribbons, which are peri- 
odically repeated structures, this develops from two par- 
tially flat bands at the Fermi level, conferring metallic 
properties [|3!|. To rationalise the above features, the 
band structure is shown in Fig. EHb), which is in good 
agreement with less approximate plane-wave DFT calcu- 
lations |63L l64| . 

The electronic states corresponding to the low energy 
bands are termed edge states. This reflects the property 
of wavefunction localisation near the edges. In fact, in 
a seminal paper of Lee and Joannopoulos |72| on the 
calculation of surface states, the honeycomb lattice with 
zigzag edges was presented as a paradigm. Other specific 
studies on tight-binding models [HJ IH3 of graphene 
revealed a single peak arising at Ep, which owes to the 



A core result of this paper relates to the fingerprints 
of the electronic band structure of the perfect one- 
dimensional crystals that correspond to the molecular 
bridges of Fig. ^andE! namely, PPE and PPP, respec- 
tively. In particular, since we study off-resonant electron 
transport within the HOMO-LUMO gap, we are not in- 
terested in the usual band structure defined for real Bloch 
vectors k. As explained below, we need instead to quan- 
tify the spatial decay of electronic evanescent states for 
energies within the energy gap E g between the top of 
the valence, E v , and the bottom of the conduction, E c , 
bands. The latter develop from the HOMO and LUMO, 
respectively, of long wires. 

When molecular junctions are formed the molecular 
electronic states hybridise with the metal wavefunctions 
giving a finite width in energy space and a renormalised 
HOMO-LUMO gap Ef (> E g ). For energies within the 
gap, matching to the metal side implies decaying wave- 
functions within the molecule. These states are the ana- 
logue of the metal-induced gap states (MIGS) in metal- 
insulator (or metal-semiconductor) junctions |55j. MIGS 
are responsible for the overall finite DOS within the Ef 
for short molecular wires because of contributions com- 
ing from both electrodes 20]. For longer wires, the MIGS 
are mainly located near the interfaces on either side of 
the molecule. Indeed, by looking at the positional de- 
pendence of the local DOS, such behaviour has been ob- 
served for Si nanowires between Al electrodes |75| and 
alkanethiol chains between Au (111) surfaces |53l . 



We come back to these points in Sec. IIVBI 

By assuming that interfacial effects at a molecular 
junction are not dominant, we turn our attention to de- 
caying states in the idealised underlying periodic systems 
of PPE and PPP. In this case, it is well known that wave- 
functions vanish exponentially in space in the forbidden 
energy domain. The decay parameter is quantified by the 
imaginary component of Bloch vectors (k = q — ik). The 
latter are obtained via the method of the complex-band 
structure, which may be understood as an extension of 
the conventional band structure to the complex plane 
[54l l55l | . Computational schemes are intimately related 
to the quest for surface states and the method is as 
follows. 

In all cases, one has to solve the Schrodinger equa- 
tion in a periodic potential, which in our tight-binding 
approximation reads 



(FT -4- H 



ika 



H 



TB 



E, S 



TB 



iTB 



ika 



iTB 



-ika\ 



-ika\ 



(10) 



Here, the matrices iff? and SfJ are defined as contain- 
ing the Hamiltonian and overlap matrix elements, respec- 
tively, between the ith and jth unit cells, a is the lattice 
constant. 

To obtain the complex-band structure, one follows the 
opposite procedure of typical band structure calculations 
in which the solution of the eigenvalue problem of the 
Hamiltonian yields the energy spectrum for any real k. 
The spectrum of k vectors associated with any real en- 
ergy E is sought. Therefore, it is useful to transform Eq. 
111)1 from an eigenvalue problem for E^ to an eigenvalue 
problem for A = e lka . This problem resembles finding 
the eigenvalues of the transfer matrix. It may be cast in 
a generalised eigenvalue problem via 



I 

\ j '(Hm,m+1 — ES mtm+ i) ^ ^ 



(11) 



which in contrast to usual formulations, can be used even 
for singular Hij — -ES;.j, as it often happens when a lo- 
calised basis set is employed. When |A| = 1, k is real; 
otherwise k is complex and we have exponentially decay- 
ing or growing solutions, which come in conjugate pairs. 

The complete band structures of perfect PPE, planar 
PPP, and non-planar PPP are calculated using Eq. ITUI 
and Eq. ^] and they are shown in Fig. EJa), 0Jb), 
and^Jc), respectively. At the right panel we present the 
conventional band structure. The spectrum of the imag- 
inary part of the complex k solutions (negative branch) 
is plotted at the left panel . Our calculated band struc- 
tures compare well with those reported for either mini- 
mal [30ll53| or plane wave [53Llr34| basis sets within other 
DFT methods in the LDA, both in the real axis and in 
the complex plane. As expected from studies of the an- 
alytic properties of the energy spectrum in the complex 
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FIG. 4: Electronic band structure for (a) PPE, (b) planar 
PPP, and (c) non-planar PPP (tetrahedral angle at 45°). The 
right panel indicates the band structure for real k in the first 
Brillouin zone, and the left panel the complex-band structure. 



plane, the complex-band structure exhibits the connec- 
tion of real bands between local extrema. Of particular 
importance for the off-resonant electron transport phe- 
nomena we discuss in the next section, is the, so-called, 
real line which joins the edges of the valence and con- 
duction bands which develop around Ep. It evidently 
corresponds to the smallest decay of wavefunctions into 
the molecular region since |k| is minimum, and it is ex- 



pected to dominate in tunnelling processes within E g . 

IV. CONDUCTANCE PROPERTIES 

A. General Features 

The local microstructure of the molecule/electrode 
interfacial contact, especially for low-dimensional elec- 
trodes [33L l3~i| . has been shown to crucially determine 
the conductance properties of a molecular junction. This 
is also seen in the transport properties of structures such 
as those of Fig. 2] To demonstrate the bare effect, we 
take the parallel and transverse geometries of a diphcnyl- 
ethynyl bridge. As shown in Fig. 03 a clear distinction 
of the conductance spectrum is evident. Compared to 
previous studies [32, |36|, |33 , a remarkable fact is that de- 
spite the carbon-carbon covalent bond at the interfaces 
and the close conformational similarity, there is a promi- 
nent dissimilarity in the conductance. The latter can be 
traced to distinct features in the local DOS for the zigzag 
and the armchair edge topologies (inset of Fig. |SJ. An 
analogue of this effect arising from the molecular bridge 
contributions to the electronic states at Ep was observed 
inRefs. [HII3- 
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E-E F (eV) 

FIG. 5: Conductance spectrum for an OPE wire with length 
N = 2. The solid line corresponds to parallel connection, 
and the dashed line corresponds to transverse connection (see 
Fig. 0. The inset shows the local DOS at the two carbon 
apex-sites for parallel (solid line) and transverse configuration 
(dashed line). 

In particular, the presence of the zigzag edge state 
has a profound influence. The conductance through the 
oligomer is non-negligible within the HOMO-LUMO gap, 
exhibiting a double resonant peak around Ep. For the 
parallel structure, on the other hand, the conductance is 
almost vanishing inside E^ 1 . The observed conductance 
resonances are in accordance with previous results show- 
ing that due to a large local DOS on the metal side for 
energies inside the HOMO-LUMO gap, an additional res- 
onant channel becomes available for through-bond con- 



duction [3(| [6(|. Indeed, the inset of Fig. [3] presents 
the local DOS at edge carbon-atoms of an isolated semi- 
infinite zigzag ribbon. The dashed line is the local DOS 
at a carbon atom (saturated by hydrogen) at the zigzag 
edge, whereas the solid line is for an atom located at 
an armchair mid-point of the terminated ribbon. The 
localised state of zigzag edge sites is featured in their lo- 
cal DOS by the characteristic spike, but is missing for 
armchair edge sites. 

To discuss generic properties of off-resonant electron 
transport properties within E g , we focus on junctions of 
the transverse connection, which exhibit additional con- 
ductance resonances within E^ due to the topology of 
the electrode/molecule interface ('worst case scenario'). 
In Fig. El we present the logarithm of the conductance 
spectrum of such carbon junctions, for various lengths 
N of the oligomers. Plot (a) corresponds to OPE, with 
valence and conduction band energies of the correspond- 
ing polymer wire E v = — 1.1 2eV and E c = 1.23eV, re- 
spectively (see Sec. IIII Bl and Fig. . Plot (b) is for 
planar OPP, with E v = -1.27eV and E c = 1.48eV, and 
plot (c) is for non-planar OPP, with E v = — 1.55eV and 
E c = 1.91eV. As it has been expected from the above dis- 
cussion, the conductance exhibits resonances around Ep, 
which are relatively high for small N. Even though g(E) 
decreases as N becomes larger, these peaks pertain for 
all wire lengths. Different oligomers, thus different elec- 
tronic configurations of the molecule and different cou- 
pling to the electrodes, do not obscure the appearance 
of these conductance resonances, indicating once again 
that they arise solely from the topology of the electrode 
contacting surface. 

The existence of both peaks in Fig. H3 can be under- 
stood by inspection of Fig. [7| There, ln(g) of an OPE 
wire of length N = 5 is plotted together with the log- 
arithm of the local DOS at the lower (dotted line) and 
upper (dashed line) electrode apex atoms in contact with 
the molecule, as a function of the energy. Each of the lo- 
cal DOS at the two contacts has a double peak inside 
the gap. The energies at which the two conductance 
peaks appear correspond to the points on the energy 
range within the gap for which there is large local DOS. 
This result can be roughly understood from the 7r-orbital 
approximation, in which the transmission is related to 
the local DOS vl/u a t the lower/upper atomic contacts 
through the relation T(E) oc v L (E)vu(E)\G 1N (E)\ 2 , 
Gin being the propa gat or between the first and the last 
sites of the molecule 1H EH S3 

From Fig. we see that the conductance decreases 
with the length of the molecule. Examining this length 
dependence in more detail, ln(g) is plotted as a function 
of N for selected energies inside the gap, in Fig. |H1 Graph 
(a) is for OPE, graph (b) is for planar OPP, and graph (c) 
is for non-planar OPP. The energy values of the curves 
defined by circle and triangle values reside in the left 
and right side, respectively, of the edge states spectrum. 
Plots indicated by the squares correspond to one of the 
two conductance resonances, and the remaining curves 
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FIG. 6: Logarithm of g versus energy E, for wires of various 
lengths TV. (a) OPE, with valence and conduction band ener- 
gies E v = — 1.12eV and E c = 1.23eV, respectively, (b) planar 
OPP, with E v = -1.27eV and E c = 1.48eV. (c) non-planar 
OPP, with E v = -1.55eV and E c = 1.91eV. 



(circles) are for an energy between the two peaks. In all 
three graphs, we evidence three types of curve. 

In Figs. IBIa) and |Hfb) and for energies away from 
those of the edge states, ln(g) falls linearly with N. This 
is reminiscent of the tunnelling behaviour, Eq. ^ In 
contrast, for energies at the vicinity of the edge state 
resonances, an anomalous length dependence appears for 



FIG. 7: The left y-axis corresponds to ln(g), and the right y- 
axis to the logarithm of the local DOS at the carbon atomic 
contacts, versus E, for OPE with N = 5. The dotted line 
is the DOS at the contact point between the lower electrode 
and OPE, and the dashed line is for the DOS at the upper 
one (see Fig. 0. 



very small N < 3, but as N increases, a linear decrease 
is also attained. We attribute the former to slight posi- 
tional shifts of the resonant energies which are evident in 
the plots as iV varies from 1 to approximately 3. After 
a complete analysis of the conductance spectrum within 
Eg , we conclude that the exponential law of Eq. ^ holds 
generally for the entire energy range within E g and for 
ultrashort molecular junctions. The observed few devia- 
tions are specific to the appearance of " unconventional" 
resonant channels induced by the electrodes. 

For the non-planar OPP (Fig. EI c ))i it appears that 
ln(g) decreases with N in an oscillating rather than a 
linear fashion. A closer look, however, reveals that one 
should distinguish between the conductance for even and 
odd N. We should bear in mind the difference in the 
interfacial contact geometry of the last benzene ring of 
the oligomer. For odd N, the benzene ring is at the same 
plane with that of the ribbon, whereas for even N it is 
positioned at an angle 9 away from the plane. This dis- 
tinction results in different coupling with the electrode, 
and therefore, in a different contact conductance as we 
analyse in Sec. IIV Bl A clear difference is also reflected 
in the structure of the conductance spectrum around Ep 
in Fig. Etc). It follows, however, that the exponential 
law applies for the two cases separately, with the same 
inverse decay length. Another viewpoint would be to 
consider an increase of the molecular length to be given 
in units of the unit cell of the non-planar PPP. 



B. Inverse Decay Length 

In the previous Section we established Eq. ^for virtu- 
ally all molecular lengths and for any energy in the regime 
of interest. Here, we identify the damping factor of the 




1.25 



FIG. 8: ln(g) versus N for various energies within the HOMO- 
LUMO gap. (a) OPE, (b) planar OPP, and (c) non-planar 
OPP. 



tunnelling conductance of an oligomer junction with the 
complex-band structure of the corresponding idealised 
polymer wire discussed in Sec. IIII Bl Consistent with 
the assumption that, to a first approximation, evanes- 
cent states in the junction would coincide with those of 
the idealised long wire and that those with the smallest 
| ac | are dominating the tunnelling process, in Fig. |5jwe 
isolate the positive branch of the real line E c — E v (solid 
line). The latter equals minus the negative branch. To- 



0.75 



0.5 



0.25 




' 

-1.25 -1 -0.75 -0.5 -0.25 0.25 0.5 0.75 1 1.25 



E-E F (eV) 



0.75 



S 

^ 0.5 



0.25 




-1.5 -1 



-0.5 0.5 

E-E F (eV) 



1 1.5 



FIG. 9: Complex-band structure (E c — E v real line) plotted 
with a solid curve and tunnelling inverse decay lengths for (a) 
PPP and OPE, (b) planar PPP and OPP, respectively. 



gether we plot the inverse decay length (3 divided by a 
factor two (dashed line), which accounts for transport co- 
efficients relating to probability amplitudes rather than 
wavefunctions. f3 is numerically obtained from a least- 
squares fit for the determination of the slope of the lines 
of Fig. |S| Graph (a) corresponds to OPE and graph (b) 
to planar OPP. For energies at which the conductance 
shows anomalous length dependence, we use only data 
for N > 2. In both graphs (a) and (b), there is a remark- 
able agreement for the entire energy window. There are 
no deviations from the relation (3 — 2 \ k{E)\ apart from 
energies at which the graphitic ribbons exhibit a large 
local DOS. 

In the discussion of Fig. 03 we mentioned that for non- 
planar OPP wires the parity of the length is an important 
element for the exact conductance spectrum, but that the 
inverse decay lengths apparently coincide for even and 
odd N. This is evident in Fig. [SJwhere the damping fac- 
tor divided by two is plotted for both cases and is com- 
pared to the complex-band structure with the smallest 
\n(E)\. The damping factors for wires with odd (dashed 
line) and even N (dotted line) are almost identical to 
each other. Furthermore, they fall on top of the posi- 
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FIG. 10: Complex-band structure (E c — E v real line) plotted 
with a solid curve and tunnelling inverse decay lengths for 
non-planar PPP and OPP, respectively. The dashed line cor- 
responds to /3/2 from OPP with odd N and the dotted one is 
for OPP with even N. 



tive branch of k which derives from the conduction and 
valence bands of non-planar PPP. Once again there is ab- 
solute quantitative agreement apart from energy points 
around singularities in the local DOS (see Fig. |3J). It is 
worth mentioning that all results for the extracted tun- 
nelling decay lengths in oligomer junctions are obtained 
for the optimised geometries, whereas, the complex-band 
structure is that of an ideal polymer wire. This shows 
the robustness of the mapping to small perturbations on 
atomic positions. 

According to the arguments of Sec lIII Bl on the wave- 
function matching between the metal and the molecular 
wire, the density of states inside the molecule should also 
reflect properties of decaying states. Namely, its projec- 
tion to an exponential should be possible, at least away 
from the electrode/wire interface. We examined this 
approach by calculating the projected density of states 
within a unit cell (MDOS) of a planar OPP wire of fixed 
length. In Fig. ITlT aL we plot its logarithm as a function 
of the monomer units i which make up N — 9, for ener- 
gies chosen inside the gap. We see that the logarithm of 
MDOS falls off almost linearly with i until roughly half 
of the molecular length is reached, and then ln(MDOS) 
increases, again in a linear fashion. Loosely speaking, the 
decreasing part of each curve corresponds to the decaying 
metal-induced gap state from the lower electrode whereas 
the increasing part comes from the gap state of the up- 
per electrode. This holds for relatively long wires as the 
one presented. For shorter molecular bridges, none of the 
contributions from each side is dominant. No apparent 
linearity as in Fig. Illf al is observed. 

To make this analysis quantitative, we calculate a 
damping factor /3mdos for the density of states, presented 
in Fig. II If . The solid line is again the smallest |k| from 
the complex band structure of planar PPP (Fig. Efb)). 
The dashed and dotted lines represent /3mdos/2 from a 
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FIG. 11: (a) Logarithm of DOS on a planar OPP wire of 
length N — 9 as a function of the monomer unit number i, 
for selected energies in the HOMO-LUMO gap. (b) Damping 
factor /3m dos for the wire of (a), from least-squares fitting of 
i = 2,3,4 (dashed line), and of i = 6, 7, 8 (dotted line). 



least-squares fit to ln(MDOS) for i — 2 to 4 and i = 6 to 
8, respectively. Both these curves agree reasonably well 
with |k|, corroborating the picture of metal-induced gap 
states. We interpret the asymmetry between the left and 
right parts of MDOS in Fig. Illf a) and (b) as differences 
on the atomic interfacial contact after optimisation. For 
example, the erroneous divergence of /3mdos/2 around 
Ep emerges at different energies in Fig. II If b') for gap 
states arising from the lower and the upper electrode, 
with the former seeing the large DOS from the upper 
electrode, and the latter seeing the large DOS from the 
lower electrode (see Fig. 0for comparison). 

We briefly comment on the correspondence between 
the inverse decay length obtained from the molecular 
density of states and the transport calculations. It is ev- 
ident that the exponential decay of wavefunctions across 
the molecular bridge as quantified by the complex-band 
structure governs both /3mdos arid /?. However, the spec- 
tral properties are determined by the diagonal Green 
function Gu . Transport is ultimately related to the prop- 
agator dj between unit cells, i.e, i ^ j. As suggested 



by the above results and the analytic properties of simple 
7r-orbital models the convergence of the latter to the long 
oligomer limit is much faster, enabling an analysis of the 
tunnelling regime via the complex k. 
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FIG. 12: Logarithm of contact conductance g D versus energy 
for all oligomer junctions under study. 



C. Contact Conductance 

The results of the previous Section elucidate that the 
inverse decay length (3 of the conductance in the tun- 
nelling regime is an intrinsic property of the oligomer 
embedded in the nanojunction. By contrast the prefactor 
go of Eq. n is largely determined by the realised interfa- 
cial contact. To investigate the so-called contact conduc- 
tance, in Fig. 1121 we present its logarithm versus energy. 
The curve denoted as OPE corresponds to the contact 
conductance of similar (or longer) molecular junctions to 
that of Fig. ^ Those curves indicated by OPP refer to 
the molecular wires drawn in Fig. [3 All data are ob- 
tained after ensuring that the exponential decay of Eq. 
n is fulfilled as explained in Sec. IIII Bl Apparently, g 
depends on the exact atomic contact (whether phenyl or 
ethynyl group to graphene) as well as the orientation of 
the phenyl rings with respect to the graphitic ribbons. 

It is possible to identify pairs of similar contact conduc- 
tances from Fig. U2 i) OPE with OPP (e), ii) OPP (a) 
with OPP (c), and iii) OPP (b) with OPP (d). Although 
the origin of the first pair is probably fortuitous, pairs ii) 
and iii) are expected from simple arguments. Moreover, 
the smaller coupling of the ethynyl-group with respect 
to the phenyl rings should give rise to the observed re- 
duced g of OPE compared to that of planar OPP. The 
different magnitude of g between non-planar oligo-para- 
phenylene for odd, OPP (c), and even, OPP (e), number 
of unit cells in the wire, expresses the importance of the 
length parity in experiments. This is in agreement with 
our previous discussion. The increased resistance for non- 
planar alignments of the phenyl end-rings with respect to 



the electrodes, e.g., OPP (b), OPP (d) and OPP (e), may 
be explained as a break-up of the conjugation. 

In addition, we estimate the contact conductance ratio 
between the parallel aligned and twisted phenyl end-ring 
geometries. We assume, as previously, only 7r-electrons 
in an orthogonal basis set. This simple model yields 
go(E) = T 2 L Tfjf(E), where, T L / V corresponds to Hamil- 
tonian matrix elements between the apex carbon atoms 
located at the wire and the graphene edges. Further- 
more, we may write (g^ /gi^) 1 ^ ~ T 45 /T° , which 
yields cos(45°) w 0.7. Considering the oversimplifica- 
tion, this agrees well with the accumulated data (inset of 
Fig. I12|) . which give a value around 0.5, apart from the 
region where the local DOS on the metal side is singular. 

Finally, we note the non-trivial energy dependence of 
the contact conductance, which may hinder simple es- 
timations based on the knowledge of j3 as those sug- 
gested in Ref. [3(|. In all our studied cases the con- 
tact conductance follows roughly the behaviour of the 
local DOS at the electrode atomic contacts, with a pro- 
nounced peak around the Fermi energy, underlining the 
importance of the electronic structure of the electrodes. 
This follows from f(E) being formally decomposed to 
vl(E)vu(E)Xm(E) in the aforementioned simple model, 
i.e., being proportional to the product of the metal local 
DOS, vl(E)vu(E) and a molecular factor Xm(E), which 
allows for resonant tunnelling at molecular orbital ener- 
gies. Although f(E) may be a slowly varying function 
of energy for leads with a seamless local DOS, especially 
deep in the tunnelling regime, it is not a priori clear that 
the off-resonant contact conductance is energy indepen- 
dent and can be determined from the on-resonance value 
[3fJ. The latter would imply equal contact conductances 
for perfectly symmetric structures based on OPP (a) and 
OPP (b) of Fig. |21 although these are at least an order of 
magnitude different (~ 0.5~ 4 ). Clearly, more detailed in- 
vestigations in diverse molecular junctions are necessary. 



V. CONCLUDING REMARKS 

We studied the generic properties of off-resonant elec- 
tron transport across molecular junctions of current 
experimental interest, formed from oligomers bridging 
graphitic electrodes. To this end, we used a compre- 
hensive analysis of the electronic structure of the various 
components in a variety of possible configurations. We 
showed how the conductance magnitude is crucially de- 
termined by the electrode spectral properties, and the 
exact topological and chemical realisation of the inter- 
facial contact. Most importantly, as a function of the 
molecular length the tunnelling conductance is almost 
always an exponentially decreasing function with the in- 
verse decay length given intrinsically by the complex- 
band structure of the corresponding idealised polymer 
chain. The properties and accuracy of the latter can be 
studied independently (cf [zij). Once the Fermi energy is 
also located with confidence, the damping factor is the- 



oretically extracted. This fact has great implications on 
comparing theoretical and experimental results. A cur- 
rent issue is the discrepancy between those. However, the 
exact contact geometry, which may have dramatic influ- 
ence, is largely unknown. Therefore, direct comparison 
of decay constants in well defined geometries, e.g., in self- 
assembled monolayers, may shed more light in our under- 
standing of electron transport across molecular junctions. 
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